Diffusion-annihilation processes in complex networks 
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We present a detailed analytical study of the A + A — > diffusion-annihilation process in complex 
networks. By means of microscopic arguments, we derive a set of rate equations for the density 
of A particles in vertices of a given degree, valid for any generic degree distribution, and which 
we solve for uncorrelated networks. For homogeneous networks (with bounded fluctuations), we 
recover the standard mean-field solution, i.e. a particle density decreasing as the inverse of time. 
For heterogeneous (scale-free networks) in the infinite network size limit, we obtain instead a density 
decreasing as a power-law, with an exponent depending on the degree distribution. We also analyze 
the role of finite size effects, showing that any finite scale-free network leads to the mean-field 
behavior, with a prefactor depending on the network size. We check our analytical predictions with 
extensive numerical simulations on homogeneous networks with Poisson degree distribution and 
scale-free networks with different degree exponents. 
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I. INTRODUCTION 

In the last few years, complex networks have become 
a new paradigm for complexity. The merging of graph 
theory together with new and classical statistical physics 
tools has lead to the development of a modern theory 
of complex network 0, , that has found fruitful appli- 
cations in domains as divers as technology (the physical 
Internet Q, the World-Wide Web 3, power grids 0), bi- 
ology (protein-protein interaction networks [fj , metabolic 
networks Q, foodwebs [3, 13)' social sciences (sexual con- 
tact networks 01 > friend ship networks , scientific col- 
laboration networks [l2lll3|| ). etc. 

The statistical analysis of many "real-world" networks 
has shown that most of these systems seem to share some 
typical features, the most relevant of them being the 
small-world property 0] and a large connectivity het- 
erogeneity, reflected in the presence of a scale-free degree 
distribution [l5j . The small- world property refers to the 
fact that, in real networks, the hop distance between two 
randomly chosen elements of the system is very small if 
compared to the total number of elements. More pre- 
cisely, if (£) is the average distance between two elements 
(or vertices), measured as the smallest number of connec- 
tions (or edges) between any vertices, and N is the sys- 
tem size (number of vertices) then, usually (£) increases 
logarithmically or slower with N. On the other hand, 
scale-free networks are characterized by a degree distri- 
bution P(k), defined as the probability that a randomly 
selected vertex is connected to k other vertices (has de- 
gree k), that decreases as a power- law, 



P(k) 



(1) 



where 7 is a characteristic degree exponent, usually in the 
range 2 < 7 < 3. For these values of the degree exponent, 
the fluctuations in the degree distribution, measured by 
the second moment (k 2 ), diverge in the infinite network 



size limit, N — > 00, giving rise a very heterogeneous con- 
nectivity lacking any characteristic degree scale. This be- 
havior is i n op position to the more classical homogeneous 
networks , which have a degree distribution decaying 
exponentially or faster and exhibit bounded degree fluc- 
tuations. 

Given that complex networks are widespread in nature, 
it becomes a quite interesting issue the characterization 
of the effects that their complex topology can have on dy- 
namical processes taking place on top of these systems. 
For example, it has been shown that heterogeneous net- 
works are remarkably weak in front of targeted attacks, 
aimed at destroying the most connected vertices Uft , 
as well as to the propagation of infective agents |19l l20j . 
These properties, which are mainly due to the critical 
interplay between topology and dynamics in heteroge- 
neous networks, are otherwise absent in their homoge- 
neous counterparts. 

Epidemic processes, chemical reactions, and many 
other dynamic processes, can all be modeled in terms of 
reaction-diffusion processes j2]| . These are dynamic sys- 
tems that involve particles of different "species" (A t , i = 
l,...n) that diffuse stochastically and interact among 
them following a fixed set of reaction rules. One's inter- 
est is usually focused on the time evolution and steady 
states of the densities of the different species pAi (t) , and 
the possible presence of critical phase transitions [22|| . 
While much is known about the behavior of reaction- 
diffusion processes on regular homogeneous lattices, the 
situation is not so well established in what respects the 
possible effects that a heterogeneous connectivity struc- 
ture can have on them. At this respect, it is interesting 
the work presented in Ref. 23], in which a numerical 
simulation analysis of the diffusion-annihilation process 
A + A — » |2^| , was performed on scale- free networks. In 
this reaction-diffusion process, particles of a single species 
A diffuse on the vertices of a network and annihilate upon 
contact (when two A particles fall on the same vertex). 
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In regular lattices of Euclidean dimension d, it is well 
known that the local density of A particles, p(x,t), is 
ruled by a Langevin equation pfij . 



dp{x,t) 

at 



DV 2 p{x, t) - 2Xp(x, tf + p(x, t)r)(x, t), (2) 



where r](x, t) is an uncorrelated Gaussian noise. Dynam- 
ical renormalization group arguments allow to show that 
the average density of A particles, p{t) — (p(x,t)), be- 
haves in the large time limit as 



p(t) Po 



(3) 



where po is the initial particle density, and the exponent 
a takes the values a = d/d c for d < d c and a = 1 for 
d > d Cl where d c — 2 is the critical dimension of this 
process. For d > d c one thus recovers the mean-field 
solution, obtained from Eq. J2J) by setting the diffusion 
coefficient D and the noise term rj(x, t) equal to zero. 

The numerical simulations of the A + A — > diffusion- 
annihilation process reported in Ref. |23l |. performed in 
scale-free networks with general degree exponent 7, gen- 
erated using the configuration model pil I27I l28l |29| . led 
the authors to conclude that the behavior in time of the 
average density of A particles can be approximated by 
Eq. Q, where the asymptotic exponent a is a decreas- 
ing function of 7 and is surprisingly larger than 1 for 
7 < 3. The authors attributed this effect to the small- 
world nature of the networks, and to the existence of hubs 
(vertices with a large number of connections). 

In spite of the potential interest of this result, no the- 
oretical arguments have been proposed so far to back up 
the numerical conclusions reached in [2j|. In this paper 
we tackle this task, by developing a mean-field analysis of 
the A + A — > process. This analysis, made in the con- 
tinuous k approximation and ins pire d in previous works 
made for epidemic spreading [T^ . l30l l3l| , results in a set 
of differential equations for the density of A particles in 
the vertices of degree k, which are valid for networks with 
arbitrary degree distribution P(k) and two vertex corre- 
lations [32|. determined by the conditional probability 
P(k'\k) that a vertex of degree k is connected to a vertex 
of degree k! [3^, 13 • The solution of these equations for 
the particular case of uncorrelated networks (in which 
the conditional probability P(k'\k) is independent of k) 
shows that, while homogeneous networks display a pure 
mean- field behavior with exponent a = 1, scale- free net- 
works with 7 < 3 in the infinite size limit exhibit instead 
an exponent depending on the properties of the network, 
i.e. a = 1/(7 — 2). Remarkably, this solution in the in- 
finite size limit shows a crossover for any finite network 
to a linear behavior 1/ p(t) ~ t, with a slope depending 
on the network size. Our analytical results are confirmed 
by means of large scale numerical simulations for both 
homogeneous and heterogeneous networks. 

We have organized the present paper as follows: In 
Sec. [n] we derive, from microscopic considerations, the 



mean-field differential equations for the A + A — > 
diffusion-annihilation process in complex random net- 
works with arbitrary degree distribution and two vertex 
correlations, quantified by means of P(k) and P(k'\k), 
respectively. We consider the case of absence of correla- 
tions, finding the density of A particles for general homo- 
geneous networks. Sec. IIIII is devoted to explicit results 
for scale-free networks, both in the infinite size limit and 
for finite size networks. In Sec. IIVI our analytic results 
are compared with extensive numerical simulations of the 
diffusion-annihilation process running on top of homoge- 
neous and heterogeneous (scale-free) networks. Finally, 
our conclusions are presented in Sec. 



II. THE A + A -> REACTION IN COMPLEX 
NETWORKS 

Let us consider the diffusion-annihilation process A + 
A — ► on a complex network of size N which is fully de- 
fined by the adjacency matrix a^-, that takes the values 
fly = 1 if vertices i and j are connected by an edge, and 
otherwise. From a statistical point of view, the net- 
work can also be characterized by its degree distribution 
P(k) and its degree correlations, given by the conditional 
probability P(k'\k). Each vertex in the network can host 
at most one A particle, and the dynamics of the process is 
defined as follows: Each particle jumps at a certain rate 
A to a randomly chosen nearest neighbor. If it is empty, 
the particle fills it, leaving the first vertex empty. If the 
nearest neighbor is occupied, the two particles annihilate, 
leaving both vertices empty. 

In order to study analytically this process in a gen- 
eral complex network, in which vertices can show large 
degree fluctuations, we are forced to consider the partial 
densities pt (t) , representing the density of A particles in 
vertices of degree k, or, in other works, the probability 
that a vertex of degree k contains an A particle at time 
t jl&l30|. From these partial densities, the total density 
of A particles is recovered from 



p(t) = J2P(k)Pk(t). 



(4) 



While it is possible to obtain a rate equation for the 
densities pk (t) by means of intuitive arguments |l9l l3l| , 
in the following we will pursue a more microscopical ap- 
proach, which can be generalized to tackle other kinds 
of problems. Let rij(i) be a dichotomous random vari- 
able taking values or 1 whenever vertex i is empty or 
occupied by an A particle, respectively. Using this formu- 
lation, the state of the system at time t is completely de- 
fined by the state vector n(t) = {ni(t), n^it), ■ ■ ■ , njv(i)}- 
Assuming that the time evolution of particles follows a 
Poisson process j2^|, the evolution of n(t) after a time 
increment dt can be expressed as 



m(t + dt) = m(t)i]{dt) + [1 - ni(t)]£(dt), 



(5) 
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where rj(dt) and are dichotomous random variables 
taking values 



If all the vertices with the same degree are statistically 
equivalent, we can set j3]| 



V(dt) = 



with probability Xdt 



1 otherwise 



a i0 n 3 {t) 



(6) 



and 



1 with probability XdtS^ Q ^ 

} k i , (7) 

otherwise 



where is the adjacency matrix and A is the jumping 
rate that, without loss of generality, we set equal to 1. 
The first term in Eq. JSJ stands for an event in which 
vertex i is occupied by a particle and, during the time 
interval (t, t + dt), it becomes empty, either because the 
particle in it decides to move to another vertex or because 
a particle in a nearest neighbor of i decides to jump to 
i, annihilating thus both particles. The second term cor- 
responds to the case in which vertex i is empty and a 
particle in a neighbor vertex of i decides to move to that 
vertex . Taking the average of Eq. (J5J , we obtain 



(m(t + dt)\n(Jt)) 



rii(t) 



n(t) - (l - 2m{t))Y^ 



dt, (8) 



equation that describes the average evolution of the sys- 
tem, conditioned to the knowledge of its state at the pre- 
vious time step. Then, after multiplying Eq. (JHJ by the 
probability to find the system at state n at time t, and 
summing for all possible configurations, we are led to 



dpi(t) 
dt 



-Pi(t) 



(9) 



where we have introduced the notation pi(t) = (rii{t)) 
and pij{t) = (ni(t)rij (t)} . 

The derivation presented so far is exact. To proceed 
further, we assume that vertices with the same degree 
are statistically equivalent That is, 

Pi(t) = pk(t) Vi e v(fc), (10) 
P ij{t) = P kk>{t) V*gV(*), j eV(fe'), (11) 

where V(k) is the set of vertices of degree k. Thus, by 
summing Eq. @ for all vertices of degree k and dividing 
by the number of vertices with this degree, 2V/., we can 
write, after some formal manipulations, 



dt 



Pk> (t) - 2pkk> (t) 1 
k' N k 



E E a * 

lEV(k) jGV(fe') 

(12) 



w k E E ^ 

iev(fc) jev(fc') 



kP{k'\k). 



(13) 



Finally, by assuming the mean-field approximation 
pkk'{t) — pk(t)pk'(t), the rate equation for the density 
Pk{t) can be written as 



dpkjt) 
dt 



-p k (t) + k[l-2p k (t)]J2 



P(k'\k) 



Pk>(t). (14) 



In the case of networks with general degree correla- 
tions, the solution of Eq. (|14fl depends on the nature of 
the conditional probability P(k'\k) and can be a rather 
demanding task |3Jj(. Therefore, in the rest of this pa- 
per we will restrict ourselves to the case of uncorrelated 
networks, in which the conditional probability takes the 
simple form P(k'\k) = k'P(k')/(k) [p. For this class of 
networks, the rate equation Eq. H14[) is simplified to the 
form 



^ = -Pk(t) + ±[l-2p k (t)]p(t), 



(15) 



where p{t) is the total density of A particles. We can ob- 
tain a differential equation for this last quantity by mul- 
tiplying Eq. I|15|) by P{k) and summing over k, namely, 



dp(t) 
dt 



-2p(i)0(t), 



where 



e W = 7^yE fcp ( fc )^W- 



(16) 



(17) 



In order to solve Eq. <|16|l . we perform a quasi-static 
approximation in the rate equation Eq. I|15|l . From the 
mean-field solution of the A + A — > process, we expect 
p(t) to be a decreasing function with a power law-like 
behavior. In this case, for large enough times, the time 
derivative of p(t) will be much smaller that the density 
proper. Extending this argument to the partial densities 
Pk(t), at large times we can neglect the left-hand-side 
term in Eq. I|15|) , and solve for pk (t) as a function of the 
density, obtaining 



Pk(t) = 



kp(t)/(k) 
l + 2kp{t)/{k)' 



(18) 



Substituting this approximation into the expression for 
Q(t), we get 



em - p{t) V fc2p(fc) 

U <fc> 2 L> l + 2kp(t)/{k)' 



(19) 
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Inserting this last expression into Eq. I|16(l . we obtain as 
a final equation for the density of A particles 



dp(t) 
dt 



= -2 



Pit? 
(k) 2 



k 2 P{k) 



2kp(t)/(k) 



(20) 



The solution of the approximate equation Eq. H2()(l will 
depend on the particular form of the degree distribution. 
The task becomes, however, quite simple for the class of 
homogeneous networks. In this case, the degree distribu- 
tion decreases so quickly that all its moments are finite. 
So, for small p(t) we can perform a Taylor expansion of 
the right-hand-side of Eq. lj2U|) , obtaining at lowest order 



dp(t) 
dt 



(k 2 



(k) 



(k)3 P[ > ' 



whose solution at large times yields 



1 

W) 



1 

Pa 



2(fc 2 ) 
(k) 2 



t, 



(21) 



(22) 



being p the initial density of A particles. This corre- 
sponds to the pure mean-field linear behavior, with a 
finite prefactor depending on the second moment of the 
degree distribution. 

Eq. I|18|l can help us to understand how this process 
will operate in heterogeneous networks. Indeed, for any 
given time t* , the partial density of vertices with degree 
larger than (k)/2p(t*) is essentially constant up to time 
t*, that is, p k (t) ~ 1/2 for k > (k)/2p(t*) and t < t* . 
The reason is that vertices with high degree (hubs) are 
more easily reached by particles than those with small de- 
gree and, with high probability they will be surrounded 
by some A particle. Then, as soon as one nearby parti- 
cle decides to move to the hub, both particles disappear 
and another nearby particle will replace the original one 
in the hub, keeping then the density constant, which is 
clearly set to 1/2. Therefore, hubs act as drains through 
which particles vanish. Besides, this process happens in 
a hierarchical way in time, that is, the partial density 
of small degree vertices decreases first, whereas the most 
connected vertices are the last to hold particles. This 
reasoning allows us to anticipate that the resulting dy- 
namics will strongly depend on the number of high degree 
vertices of the network and, consequently, on the degree 
distribution. 



III. SCALE-FREE NETWORKS 

For heterogenous networks with a diverging second mo- 
ment, as in the case of scale-free networks, we have to 
consider more carefully the solution of Eq. Let us 

focus in generalized uncorrelated scale-free networks in 
the infinite size limit, which are completely determined 
by the normalized degree distribution 



P(k) = (7- \)m'<- 1 k-' 1 , 



(23) 



where 2 < 7 < 3 is the degree exponent, m is the min- 
imum degree in the network, and we are approximating 
k as a continuous variable. The average degree is thus 
(k) = 777,(7 — l)/(7 — 2). For a scale-free network in the 
limit N — > 00, Eq. (|20|l can be written as 



dp(t) 
dt 



-2p(t) 



2(7 



l)m 



7-1 



p(<)F(l,7-2,7 



k 2 ^ 



1 , 2fc P (t; 

(7-1) 
2( 7 ~2)p(t) 



-dk 



(24) 



where F is the Gauss hypergeometric function [3_^ 
Therefore, the density as function of time can be im 
plicitly expressed in terms of the integral 



t = 



p(ty 



d,z 



I'n 



zF(l, 7 -2,7-1, - 



' 2(-y-2) I 



(25) 



For very large times and small densities, we can use the 
asymptotic expansion of the Gauss hypergeometric func- 



tion, F{1, 7 — 2,7—1, — z) ~ z 



2-7 



00, to obtain the 



scaling behavior of the density with time, 



P(t) 



t 



1/(7-2) 



(26) 



This same result can be derived in a more intuitive fash- 
ion starting from the argument presented at the end of 
the previous Section. If hubs act as drains of particles, 
then the rate of change in the total density of particles 
will be proportional to the density of hubs. More pre- 
cisely, identifying the relevant hubs at time t as those 
vertices with degree larger than (k)/2p(t), we have that 



dp(t) 
dt 



P(k)dk ~ p(t) 



7-1 



(27) 



from where we obtain the same trend given by Eq. 12611 . 
As 7 decreases the probability of having high degree ver- 
tices as drains in the network increases. From Eq. I|26|) 
we see as well that when this happens, the decrease in 
the density gets faster. Therefore, the bigger is the de- 
gree of the hubs in the network, the faster is the process 
of annihilation. 

We can obtain an exact solution for p(t) in scale-free 
networks with degree exponent 7 = 3. In this case, we 
have P(k) — 2m 2 k~ 3 and (fc) = 2m. The differential 
equation for the particle density reads now 



dp(t) 
dt 



-P(tf 



1 + kp(t)/m 



dk 



(28) 



whose solution is 

[M- 1 dz 
l p -i ln(l + z) 



t 



L i (l + p(t)^ 1 )-L i (l+ Po - 1 ), 



(29) 
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where Li is the logarithmic integral function |36|. For 
small p(t), we can exploit the asymptotic expansion of 
the logarithmic integral function, L;(z) ~ zj lnz, z — > oo, 
to obtain 



1 



tint. 



(30) 



That is, in the infinite network size limit, the particle 
density of a scale-free network with 7 = 3 follows the 
same decay as the mean-field solution, with a logarithmic 
correction. 

The conclusion of this analysis is that the particle den- 
sity in the A + A — ► process in a infinite size scale-free 
network with degree exponent 7 follows the asymptotic 
form at large times 



1 



P(t) 



with characteristic exponents 
0(7) = 



1/(7-2) 
1 



2 < 7 < 3 
7 > 3 





0(7) = < 1 




2 < 7 < 3 
7 = 3 
7 > 3 



(31) 

(32) 
(33) 



The results for 7 > 3 are a natural consequence of the 
lack of degree fluctuations ((fc 2 ) < 00) in this kind of 
scale-free networks. 

The analytical exponents determined above, however, 
are strongly affected by finite size effects. Indeed, for a 
power-law degree distribution, the largest weight in the 
sum in Eq. H20(l is carried by the large fc values. If the 
network is composed by a finite number of vertices, N, as 
it always happens in numerical simulations, it has a cut- 
off or maximum degree k c (N), which is usually a function 
of the network size |37j . Thus, there exists a cross-over 
time t c , defined by 



2k c (N)p(t c ) 
(k) 



« 1, 



(34) 



such that, for t > t c the particle density is so small that 
we can approximate 



e(t) 



(fc) 2 



Y,k 2 P(k)=p(t) 



(35) 



In this time regime, we will observe an effective linear 
behavior for the inverse particle density, 



1 



2(fc 2 ) 
(fc) 2 



(36) 



whose effective slope depends of the network size as 

k c (N) 

(fc 2 ) = V k 2 P(k) ~ fc c (A0 3 " 7 . (37) 



That is, in finite size scale-free networks, the density of 
A particles at very large times has the same behavior as 
in homogeneous networks, i.e. \/p{t) is linear in t, but 
now with a slope that depends on the network size N. 
Thus, the exponents larger than one, given by Eq. 13211 . 
can only be numerically observed as transients in very 
large networks. 



IV. NUMERICAL SIMULATIONS 

As a check for the analytic predictions developed in the 
previous Sections, we have performed extensive numer- 
ical simulations of the diffusion-annihilation process on 
top of different model networks, both homogeneous and 
heterogeneous. The simulations are implemented using a 
sequential updating scheme j22j . An initial fraction po of 
vertices in the networks is randomly chosen and occupied 
by an A particle. At time t in the simulation, a vertex 
is randomly chosen among the n vertices that host an 
A particle at that time. One of its neighbors is selected 
also at random. If it is empty, the particle moves and 
occupies it. If it contains a particle, this is annihilated 
with the one in the first vertex, becoming both vertices 
empty, and the number of occupied vertices is updated 
n — > n — 2. In both cases, time is updated as t — > t + 1/n, 
where n is the total number of particles at the beginning 
of the simulation step. In all the simulations presented in 
this Section, we set the initial particle density po = 0.5. 



A. Homogeneous networks 

The class of homogeneous networks refers to networks 
with a degree distribution peaked at an average value 
(fc) and decaying exponentially or faster for fc 3> (fc) and 
fc <C (fc). Typical examples of such networks are the 
Erdos and Renyi model [33 and the small-world model 
proposed by Watts and Strogatz (WS) [Tj|- We will focus 
in the latter in order to perform computer simulations. 
The WS model is defined as follows [HE^: The start- 
ing point is a ring with N vertices, in which every vertex 
is symmetrically connected to its 2m nearest neighbors. 
Then, for every vertex, each edge connected to a clock- 
wise neighbor is rewired to a randomly chosen vertex with 
probability p, and preserved with probability 1 — p. This 
procedure generates a graph with a degree distribution 
that decays faster than exponentially for large fc, and av- 
erage degree (fc) = 2m. We will consider the WS model 
with p — 1, that is, in which all edges have been rewired. 
In this limit, the degree distribution of the WS network 
takes the form 391 



P(fc) 



,fe— m 



(k — m)\ 



(38) 



Therefore, all its moments (fc n ) are finite for any value 
of n, and we should expect Eq. IL'l'l) to provide a good 
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FIG. 1: Inverse average particle density p(t) as a function of 
time for the A + A — > process in WS networks. The dashed 
line corresponds to a behavior l/p(t) ~ t. 



A 



process 



approximation for the dynamics of the A- 
for p(t) 0. 

In order to check this fact, we have carried out large 
scale numerical simulations of the A + A — » reaction on 
WS networks with p = 1 and minimum degree m = 2, 
which corresponds to an average degree (fc) = 6. Simu- 
lations were performed on graphs of size N = 10 6 , aver- 
aging over 1000 reaction processes over 10 different real- 
izations of the random network. In Fig. ^ we show the 
total density of A particles as a function of time. As 
we can observe, the inverse particle density follows quite 
precisely a linear behavior with time. A linear fit for the 
last two decades of t values yields a slope ps 1.50, which 
is not too different from the predicted value in Eq. ( 1221 . 
2{k 2 )/(k) 2 = 2.25. 



Heterogeneous networks: The Barabasi-Albert 
model 



The Barabasi-Albert (BA) graph was introduced as 
the first growing network model capable of producing 
as an emerging property a power-law degree distribution 
[l5| . This model is based on the preferential attachment 
paradigm, a rather intuitive mechanism in which new 
individuals tend to develop more easily connections with 
already well connected individuals. The model is defined 
as follows: We start from a small number m,Q of vertices, 
and at each time step, a new vertex is introduced, with m 
edges that are connected to old vertices % with probability 



n(fei) 



(39) 



where ki is the degree of the ith vertex. After iterat- 
ing this procedure a large number of times, we obtain 
a network composed by N vertices, minimum degree m, 
and average degree (k) = 2m, with a degree distribution 



FIG. 2: Check of Eq. in BA networks. The networks 
have size N = 10 5 . 




FIG. 3: Inverse average particle density p(t) as a function of 
time for the A + A — » process in BA networks. The full line 
is the analytic solution Eq. 11291 : the dashed line corresponds 
to the finite size behavior 1/ pit) ~ t. 



(k) = 2m 3 fc -3 and almost vanishing degree correlations 
In the results presented here we consider the param- 
eters uiq = 5 and m = 2, corresponding to an average 
degree (k) = 4. Simulations were performed on networks 
of size N — 10 5 and N — 10 6 , averaging over 1000 runs 
in 10 different network samples. 

A first issue that we have explored in this network is 
the validity of the quasi-stationary approximation made 
to obtain Eq. (|18f) . If this approximation is valid, we 
should observe that the quantity 



k 



1 



Pk(t) 



2 = 



{kl 
P(t) 



(40) 



is independent of k. We have performed this check in 
Fig. |21 As we can observe, the curves collapse quite 
nicely for t larger that to ~ 10. Therefore, we should 
expect that the results obtained from the approximation 
in Eq. (|20fl will be valid above this characteristic time. 

In Fig. |3| we represent the inverse particle density as 
function of time. In this plot we observe that, except 
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for very early times, the initial time dependence is very 
accurately described by the analytic solution Eq. 
which can be approximately described by a linear be- 
havior with a logarithmic correction, l/p(t) ~ tint (full 
line). At larger times, the function l/p(t) crosses over 
to a linear dependence (dashed line), showing the onset 
of finite size effects. The poor agreement at early times 
should be attributed to the approximation in Eq. <|20|) . 
which is valid only for times larger than to ~ 10. 



C. Heterogeneous networks: The uncorrelated 
configuration model 

The construction of uncorrelated scale-free networks 
with arbitrary degree exponent is a nontrivial issue. Tra- 
ditionally, this kind of networks were gen erated using the 

configuration model (CM) HHHm^, in which ' start " 
ing from a degree sequence extracted according to the 
desired degree distribution, edges are randomly created 
between vertices, respecting the preassigned degrees. It 
has been shown, however, that for scale-free distributions 
with diverging second moment, and in the absence of 
self-connections (a vertex joined to itself) and multiple 
connections (two vertices connected by more than one 
edge), this procedure generates in fact degree correla- 
tions H3, EI- The origin of this phenomenon can be 
traced back to the effects of the cut-off, or maximum ex- 
pected degree, k c (N), which must scale at most as A 1 / 2 
in order to allow for a full random uncorrelated edge as- 
signment |4^ . 

To solve this question, it has been recently proposed 
the uncorrelated configuration model (UCM) |43j, that is 
defined as follows: 

1. Assign to each vertex i in a set of A initially dis- 
connected vertices a degree ki, extracted from the 
probability distribution P(k) ~ /c~ 7 , and subject 
to the constraints m < ki < A 1 / 2 and ^T. ki even. 

2. Construct the network by randomly connecting the 
vertices with edges, respecting the pre- 
assigned degrees and avoiding multiple and self- 
connections. 

Using this algorithm, it is possible to create scale-free 
networks whose cut-off scales as k c (N) ~ A 1 / 2 for any 
degree exponent 7, and which are completely uncorre- 
lated. The results presented below were obtained from 
simulations performed on networks generated from the 
UCM algorithm with minimum degree m = 2 and sizes 
ranging from A = 10 5 to A = 10 6 . Averages were per- 
formed over 1000 runs in 10 different network. 

In Fig. 2| we plot the inverse particle density from com- 
puter simulations in networks with different degree expo- 
nent 7. From this plot we observe that, at the initial time 
regime, the growth of this function is faster for smaller 
values of 7, in agreement with the theoretical prediction 
Eq. (|2(j|l . At larger times, on the other hand, finite size 
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FIG. 4: Inverse average particle density p(t) as a function 
of time for the A + A — > process in uncorrelated scale-free 
networks with different degree exponents. The dashed line 
corresponds to the finite size behavior l/p(t) ~ t. 

affects take over, and we observe again a linear regime 
(compare the slope with the reference dashed line in the 
Figure). For 7 < 3, however, we cannot compare di- 
rectly with the full analytic prediction Eq. (|2ul) . as we 
did for the BA model. The reason is the following: From 
Eq. (|34f> . we have that the density at the crossover time 
scales as 

p(t c ) ~ k c {N)-\ (41) 

Then, from Eq. I|26|) . we can estimate 

t c ~ k c {Ny~ 2 . (42) 

That is, the crossover time t c diminishes with decreasing 
7. It turns out that, even for 7 = 2.5, t c is so small that 
we cannot observe the transient regime l/p(t) ~ i 1 ^ 7-2 ), 
and only the finite size effects, with its associated linear 
behavior, are apparent. System sizes larger than those 
available for our present computer resources should make 
visible this crossover phenomena. 

Finally, we have checked the slope dependence of the 
linear behavior given by the finite size effects. Accord- 
ing to Eqs. (|3"rJjl and (|oT|) . we should observe, for a fixed 
degree exponent, a slope that grows with the system size 
as /c c (A) 3 ~ 7 ~ jv( 3 ~ 7 )/ 2 , where we have used the scal- 
ing of the cut-off given by the UCM algorithm. We have 
analyzed this point in Fig. [5] For a degree exponent 
7 = 2.5, we observe that the curves for increasing values 
of A show a correspondingly increasing of the slope in 
the final linear region. We have estimated this slope by 
performing a linear fitting to the whole curve (inset in 
Fig. 0. This slope increases approximately as A 1 / 4 , in 
very good agreement with our theoretical predictions. 

It is noteworthy that these finite size effects, so re- 
markable in the simulations presented in this Section, are 
however not evident in the numerical work developed in 
Ref. [2^ ■ The reason of this fact could be that the simu- 
lations in |23| were performed in networks generated with 
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FIG. 5: Inverse average particle density p(t) as a function 
of time for the A + A —* process in uncorrelated scale-free 
networks with degree exponent 7 = 2.5 and different network 
sizes. Inset: slope in the linear regime as a function of TV 1 / 4 . 
The good linear behavior validates the finite size solution in 
Eq. fljl. 

the CM algorithm. These networks, as we have discussed 
above, exhibit some degree of correlations 0, 0, Elij . 
and therefore could be outside of the regime of validity 
of the analytical calculations we have presented here. 

V. CONCLUSIONS 

In this paper we have presented a detailed analytical 
study of the A + A — > diffusion-annihilation process in 
complex networks 23]. For uncorrelated homogeneous 
networks with bounded degree fluctuations, we recover a 
behavior compatible with the standard mean-field solu- 
tion of the process, that is, the inverse particle density 
grows linearly with time, l/p(t) ~ t, with a constant 
prefactor. In the case of uncorrelated heterogeneous net- 
works, characterized by a scale-free degree distribution, 
we observe instead that, in the infinite network size limit 



TV — ► 00, the inverse particle density increases as a power 
law with time, \/p(t) ~ t a ^\ with an exponent larger 
than one and that depends on the level of heterogeneity 
of the network, i.e. 0(7) = 1/(7 — 2) for 7 < 3. For 
7 = 3, a linear growth with logarithmic corrections sets 
in, while for 7 > 3 we recover again the mean-field so- 
lution typical of homogeneous networks. In the case of 
scale-free networks, we have also analyzed the effects of 
a finite network size in this dynamical process. For any 
value of 7 < 3, we have shown that, at large times, the 
inverse particle density in a finite network crosses over to 
a linear behavior, with a slope that is an increasing func- 
tion of the network size TV, i.e. l/p(t) ~ N&'^/H. In 
order to check our results, we have performed extensive 
numerical simulations on both homogeneous and hetero- 
geneous networks, obtaining a convincing evidence for 
our predictions. 

An interesting conclusion that can be extracted of the 
present work is the very strong incidence that finite size 
effects can have on dynamical systems on top of scale- 
free networks. While these size effects have been already 
discussed in the context of epidemic spreading 0] j the 
v4+j4 — > process analyzed here provides a further exam- 
ple, in which analytical results and numerical simulations 
show a striking agreement. 
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